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A8STRACT 


An Advanced Simulation of Digital Filters has been 
implemented on the !BM 360/67 computer utilizing Tektronix 
hardware and software. The program package is appropriate 
for use by persons beginning their study of digital signal 
processing or for filter analysis. The ASDF programs provide 
the user with an interactive method by which filter pole and 
zero locations can be manipulated. Grpahical output on both 
the Tektronix graphics screen and the Versatec plotter are 


provided to observe the effects of pole-zero movement. 
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lf. $tNTRONUCTION 


Digital methods have recently begun to dominate’ the 
selection of transmission and filtering schemes. When the 
student first confronts the discrete domain, he encounters a 
series of unfamiliar concepts. While most of the discrete 
concepts have parallels in the continuous or s-domain, this 
dualism [fs not apparent initially. 

The programs comprising the Advanced Simulation of 
Digital Filters (ASDF) provide the user with a convenient 
method for developing intuition jin the discrete domain. 
This is achieved by allowing the user to generate digital 
filters by directly entering the ooles and zeros into a 
graphics terminal. The ASDF programs then perform. the 
calculations required to display the characteristics of the 
user generated filter. 

The overriding objective of the ASDF program package is 
to provide the user with an interactive method for modifying 
the filter pole/zero locations in the z-plane, as well as to 
compute the filter responses. Both time and frequency 
responses are computed and displayed graphically on the 


Tektronix 4912 graphics terminal. [Ref. 9] 





Pe ce eeeciVOn Or ASDE ALGORITHMS 


There are several parameters which are particularly 
important in selecting the algorithms used in the ASDF 
programs. Since one of the primary objectives of these 
programs is to provide the user with clear explanations of 
methods for computing digital filter characteristics, direct 
or straightforward methods should be used. The importance 
of developing a user's intuition for the discrete domain and 
digital filters Suggests that graphical rather than 
numerical outputs would be most appropriate. The use of 
graphical outputs to characterize a filter's responses 
implies minimum importance of precision in these results. 
Since the reason for computing the filter characteristics is 
to develop the user's intuition for the discrete domain, 
there are no pressing constraints on the selected algorithms 
with regard to computational efficiency, speed, or memory 
space utilization. In fact, the final ASDF ovackaze spends 
the greatest portion of its execution time making the 
interface with the z-plane, and the form of the output easy 
to understand and to use. The greatest part of the memory 
space utilized for the ASDF programs is consumed in_~ the 
interfacing processes rather than in’ the methods” by which 


the data is stored, used, or computed. 





PeeeeeALGORITHMS SELECTED 

Because of the above stated objectives, a simple 
algorithm was selected for computation of time domain 
responses of digital filters. Once the user enters’ the 
desired pole and zero locations, a direct-form-two 
difference equation representation of the filter ins 
generated [Ref. 5]. The unit sample sequence is’ used for 
the filter input, the difference equations are solved 
recursively in the discrete time domain and the unit sample 
response is recorded as the output. The unit step response 
is calculated by a similar method. 

It is known that the direct~form-two filter 
implementation is very sensitive to parameter variations, 
and as such cannot be utilized in situations requiring 
extreme accuracy. The advantages associated with simplicity 
of understanding the algorithms, the ease of programming the 
general case, and the deemphasis of accuracy in outputs for 


marginally stable filters outweighs this disadvantage. 


B. DEFINING THE DATA STRUCTURE 

A detailed description of the important variables in the 
ASDF programs is given in Apnendix 8B. There is, however, an 
overiding data structure common to al}! portions of the ASDF 
package which is important. 

The array POLZRO holds the locations of the poles and 


zeros of the filter under analysis. In all cases, zeros are 





stored in POLZRO(I,J) for values of I from one to ten. Pole 
ieeations are stored in POLZROCI,J) for values of I from 
eleven to twenty. Only the top half plane pole of = each 
complex pair is stored. There are five sections of 
POLZROC( I,J), i.e. values of J from one to five. POLZROCI,1) 
is the real part of. the rectangular representation of the 
en root. POLZRO(!1,2) is the imaginary part of the 
rectangular representation of the ith root. POLZROCI,3) and 
POLZRO( 1,4) are the magnitude and angle portions of the ith 
root location in polar form. POLZRO(!I,5) is the root tyne 
designator. Root type designators are: one for real zeros, 
two for complex zeros, three for real poles, and four for 
complex poles. 

The array I!ROOTS(I) holds the numbers of roots of each 
type of root already entered into the system. IROOTS(1) and 
me90O1S(3) hold the number of real zeros and poles, 
respectively. JIROOTS(2) and IROOTS(4) hold the number of 
complex zero and pole pairs, respectively. 

The variable IERROR is the error flaz, which is” set 
whenever the user generates an error. There are numerous 
possible errors. These include trying to delete roots which 
have not been entered, or adding roots which cause the 


maximum system order to be exceeded. 
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Pou tes IMPLEMENTED SOLUTION 

Having determined the general form of the solution it 
remained to generate the software implementation. There 
were five major areas to be programmed: the interface with 
the z-plane, the time response, the frequency response, the 


plotting of output curves, and the supervisory routines. 


Peer NTERFACE WITH THE Z-PLANE 

Since the method by which the user modifies the pole and 
zero locations in the z-plane is of no theoretical interest, 
the discussion of this implementation will be brief. 
Extensive coverage of the subroutines employed and_ the 
definitions of the variables used is reserved for Appendix 
5. 

The ASDF command which allows the user to make pole-zero 
modifications in the 2z-plane is POLZRO. This command 
initiates execution of a load module which is structured as 
depicted in figure 3-1. Each of the boxes shown represents a 
subroutine. Lines with arrow heads on both ends indicate 
that program control passes to the subroutine and_ then 
directly back to the calling routine with no intermediate 
branches. The main orogram, MAIN, initializes the required 
arrays, and tests to ensure that the user has’ generated a 


realizable causal filter. 


shat 





ow LP 


MAIN 
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Figure 3-1 Flowchart Of POLZRO 





The actual interaction with the z-plane is accomplished 
in subroutine UNTCIR. The structure of UNTCIR is diagrammed 
in figure 3-2. The large double ended arrow indicates that 
there are many subroutines of the form CRRMMM, The 
subroutines LETTER and BOXUC generate the unit circle and 
command array. PLOT-10 software CRef. 10] statements 
determine the location of the user positioned cursor. This 
cursor position is adjusted with software tabs to one of 
twenty-six possible locations. Based on the X and Y 
location of the adjusted cursor position a specific 
subroutine is selected for execution. The Flashing 
alphanumeric cursor is positioned in the selected command 
box using the adjusted cursor location as a reference. The 
subroutine then reads two alphanumeric characters typed by 
the user. The combination "space", "e'" indicates that the 
user wishes the selected subroutine to be executed. The 
combination "space", "x" indicates that the user wishes to 
terminate the pole zero manipulation phase of the program. 

All of the commands which change the status of the 
z-plane have names of the form CRRMMM where C is either "A" 
for add root(s) or "D" for delete root(s). The RR indicates 
root type: RZ - real zero, CZ - complex zero, RP - real 
pole, CP -=- complex pole. The MMM stands for the mode of 
data entry. 1AG stands for interactive graphics i.e. the 
roots will be adjusted with the cursor. POL and RCT stand 


for polar and rectangular, respectively. The general 


les 





Iie a EAR 


OOF bik BOXUC 


ERROR UNTCIR I EXOMMUKG 


CRM! 


Figure 3-2 Structure Of Subroutine UNTCIR 
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structure of the CRRMMM subroutines is described in figure 
Bea « 

All CRRMMM subroutines follow the same general pattern. 
The POLZRO table is checked to determine whether or not the 
issued command is valid. Poles and zeros must be present 
before they can be deleted. There cannot be more than ten 
poles or ten zeros. The desired pole or Zero location is 
next acquired. For subroutines adding roots! the entered 
location is. that of the root(s) to be added. For 
subroutines deleting roots the entered location is any 
location near the root(s) to be deleted. Using subroutine 
LOCATE, the POLZRO table is searched for the root of the 
correct type which ts nearest the snecified location. This 
nearest root is then highlighted. In all cases, except when 
adding roots in the !AG mode, the user is asked whether or 
not the correct location has been determined. Once the root 
location is verified the root is added and nlotted with 
subroutines STRRTS and PLOTRT. For subroutines deleting 
roots the POLZRO table is adjusted with subroutine UPDATE 
and the unit circle is redrawn. 

The user may make anerror in carrying out’ the above 
Steps, t.e., by placing a pole outside the unit circle, by 
exceeding the maximum system order, or by entering a zero 
outside ten units from the origin. In each of these cases an 


1 ereafter, poles and zeros will be refered to as roots. 
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Figure 3-3a Structure of CRRMMM Subroutines 
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Figure 3-3b Structure Of CRRMMM Subroutines Continued 
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Figure 3-3c Structure Of CRRMMM Subroutines Continued 
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error flag is set. Subroutine ERROR handles the printing of 
messages and subroutine LST prints a listing of the current 


filter status before the program resumes. 


B. TIME RESPONSE 

The transfer function of a digital filter is defined in 
equation 3.1. The Ze Dae Aye and b. will be defined ina 
subsequent section. The variable A is an arbitrary constant 


assigned by the user and is called the filter gain. 


M M N 
fay (z-z.) eke 
j=] - pe0e 
H(z) = ———_—____———_- > ———__—__—_—_ (3,1) 
N N k 
TT (z-p.) joss re 
i=] = ksi * 


[Inherent in generating the time domain response with a 
direct-form-two filter implementation are the determination 
of the coefficients ay and Die im equation (3-1), and the 
computation of the unit sample and unit step sequences. 
Having determined the filter coeficients, the time domain 

response of the filter is determined by using the 
appropriate filter input sequence and recording the output 
of the filter for the desired number of points. This is 
accomplished by recursively solving a PCR erence equation 
which characterizes the filter for 1924 time points. The 


standard direct-form-two filter implementation is shown in 


figure 3-4. The 2 terms represent unit time delays. This 
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figure corresponds to the difference equations: 


N N 
y(n) = f ay (n-k) + 5 
k=] r=0 


bx Cn-Ner) (362) 

To compute the filter's time response, three operations 
are required: calculate the ay and 2 coefficients, generate 
the input sequence x(k), k=0, 1, . . . 1024, solves the 
difference equations and stores the filter output values. 
There are four central subroutines which perform these 
operations: ASUBK, BSUBR,  RSPNSE, and COEFF. Subroutines 
ABSUK, BSUBR, and COEFF generate the filter coefficients 
starting with the pole and zero locations specified. 
Subroutine RSPNSE solves the difference equations by 
coordinating the other three routines, generating the 
appropriate input sequence, and storing the output sequence 
values. 

The ASUBK and _  8SUBR subroutines are identical in 
structure and purpose. ASUBK computes the degree of the 
denominator polynomial as two times the number of complex 
pole pairs plus the number of real poles. A comnlex vector 
is formed named POLES, which holds a listing of all the 
poles of the filter under analysis. The POLZRO array stores 
real poles in the form X+#j30. Only the top half plane pole 
of each complex pole pair is stored, i.e., (X+jY¥), where X 
and Y are non-negative. The vector POLES has’ entries for 


all real poles and both poles of each complex opole pair. 
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Once the POLES vector has been generated the subroutine 
COEFF multiplies out the factors of the form (z-p.) to 
generate the coefficients of the denominator polynomial. 

The algorithm for multiplying out the denominator 
polynomial ts essentially the same as would be done by hand. 
multiplication of aseries of Ilfiteral factors of the form 
(2-p.). The vector POLES fs passed to the COEFF subroutine 
as the vector V. Since the denominator factors are of the 
form (z-p.), the V vector is changed to hold the values “Pp, 
lf the system is of order less than ten, the remaining 
positions of V are filled with zeros. The degree of the 
denominator polynomial [ts subsequently used to determine how 
many poles are actually located at the origin. 

The vectors used in the subroutine COEFF hold the 
coefficients for the polynomial heing senerated. The 
literal powers of "z"' are supressed. The vector H1 always 
contains the partial product of the K factors oreviously 
multiplied together. The vector H2 represents the old 
waeetal product times the "z" of the (z-p;) factor. Since 
the coefficient of "z'' is always one, H2 is a left shifted 
version of H1. The vector H%3 represents the partial product 
times the next pole location (p;). This mext -p; value is 
stored in V(1). The new H1 vector is formed by adding the 
H2 and H3 vectors. H1 now contains the product of the first 
fea factors. The values in the VY vector are next shifted 


such that the next pole location to be multiplied is always 
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in the V(1) postition. These four steps are repeated ten 


times to menerate the coefficients of the complete 


polynomial in Z The formula for the transfer function 
requires coefficients of the denominator polynomial in terms 
of z~+ so the order of the coefficients is reversed. 

Figure 3-5 depicts several iterations of the algorithm 
used in the COEFF subroutine. Figure 3-5a shows the initial 
conditions for each of the vectors defined in the above 
paragraph. BicvgemooscmeshOwse the Svyalites of the vectors 
after the first iteration. Figure 3-5f shows the vector 
values after the second iteration. The third and subsequent 
iterations follow this pattern. 

In precisely the same way as7 discribed above, the ay. 
coefficients are senerated with the ASUBK subroutine. The 
zero locations are stored in POLZRO exactly as are the 
poles. The generation of the vector ZEROS and the use of 
the subroutine COEFF parallels the BSUBR subroutine exactly. 

The significant variables in subroutine RSPNSE are: 
XSUBN, INPUT, OUTPUT, FLTRGN, A, and 8. The vectors A and 8 
hold the already computed values of a, and b,. AC 1) holds 
@gz AC2) holds a, . . . AC11) holds ajg. The b, are stored 
similarly in the B vector. XSUBN is the current value of the 
input sequence. INPUT is the value of the node so labeled 
in figure 3-4. OUTPUT is arly noted on figure 3-4, 


FLTRGN fs the value of the filter gain, and is assigned by 


the user. The vector Z holds the delayed INPUT values. 
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Figure 3-5Sa. Initial Conditions Of COEFF Vectors. 
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Figure 3-5b. Values Of COEFF Vectors Before The End Of The 


First Iteration. 
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Figure 3-5c. Values Of COEFF Vectors starting The Second 


fteration. 
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Figure 3-5d. Values Of COEFF Vectors before The End Of The 


Second Iteration. 
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Figure 3-5e. Values Of COEFF Vectors Starting The 


Third Steration. 
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Figure 3-5f. Values Of COEFF Vectors before The End Of The 


Third ![teration. 
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mmeeis no delay; Z(2) is delayed one time unit; 2Z(3) jis 
delayed two time units, etc. 

The subroutine RSPNSE generates the appropriate input 
sequence. For the unit sample response the input sequence 
mee 0, 0, 0, ... . >%FFor the unit step response the input 
sequence is 1, 1, 1, l, ... . The appropriate value of 
XSUBN is determined by the tyne of input sequence. The 
variable INPUT is then computed to be the sum of XSUBN and 
the a's times the delayed INPUT values. The filter output 
is the sum of the bs times the delayed INPUT values. 
These output values are then multiplied by the filter gain, 


FLTRGN, and stored in the appropriate array: YI for the 


unit sample response, YS for the unit step response. 


Gy FREQUENCY RESPONSE 

For discrete systems, the frequency response is found by 
evaluating the system transfer function at points on the 
unit circle. Subroutine FROANCY generates both the ohase 
and the magnitude of the transfer function. 

The important variables in the FROQNCY subroutine are: 
AA, BB, THETA, Z, NNUM, and NDEN. The variables AA and BB 
are the double precision values of the A and B- vectors 
described in the section on time response. THETA is the 
angular displacement around the unit circle for which a 


value of the frequency response is calculated. Z is the 
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classical z-plane variable: 
Z=exp(j THETA) C325) 


NNUM is the number of numerator terms. NDEN is the number 
of denominator terms. 

The frequency response is computed for 1024 points 
around the unit circle. All points computed are for THETA 
greater than zero and less than wT. The values of Z for 
THETA greater than -7 and less than zero are neither 
computed nor plotted since the magnitude of the transfer 
function is an even function of THETA and the phase response 
is an odd function of THETA. 

The generation of the ay. and b. coefficients has already 
been discussed in the section on time domain response. The 
evaluation of H(z) is straightforward. The value of THETA is 
incremented through 1924 steps from zero to. For each 
value of THETA generated, the value of Z is computed. The 
value of H(z) is then computed in two parts. The numerator 
Peeene sum of NNUM terms of the form 5b, times 2” for 
r=0, 1, 2,... . NNUM-1°. The denominator is one minus the 
sum of (NDEN-1) terms of the form ay times Das for 
mere 2, 3, 2 © NDEN-1°. Dividing numerator by denominator 

al 


NNUM is equal to the degree of the numerator plus one. 


2 NDEN equals the degree of the denominator plus one. 
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yields H(z). The phase of H(z) Is then determined using the 
principal values of the arctangent function; the results are 
stored in the vector YP. The magnitude of H(z) is computed, 
multiplied by FLTRGN, the filter gefain, and subsequently 


stored in the vector YM. 


Pee OUTPUT PLOTS 

The output plots generated by ASDF are generated by two 
different software programs each being different 
implementations of the same plotting requirements. One set 
of plots is generated using PLOT-10 software on_ the 
Tektronix 4012 terminal. Another set of plots utilizes the 
Versatec plotter and associated software. [Ref. 7] The 
PLOT-10 software is adaquately covered in reference 10. The 
Versatec software is explained in Technical Note-034 in the 
computer center. 

The plots which are generated on both devices’ are 
essentially the same. There are five plots generated on the 
4012 terminal: unit sample response, unit step response, 
phase response, magnitude of the transfer function, and 
magnitude of the transfer function in decibels. Both of the 
time domain plots are terminated when the output value has 
stabilized to within one percent of the maximum generated 
system output. All 1024 points of the phase and magnitude 


plots are plotted. 
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The form of the Versatec output is similar to that of 
the 4012. One additional plot is included. It shows the unit 
circle with poles and zeros of the filter being analyzed. 
Plotting is terminated when the filter output stabilizes to 
within one percent of the maximum system output generated. 

On both devices, the values of unit sample response are 
stored in a vector YI and plotted with subroutine PLTIMP. 
The unit step response values are stored in YS, and plotted 
with subroutine PLTSTP. The phase response is stored in 
vector YP, The magnitude of the transfer function is stored 
in vector YM. The magnitude of the transfer function’ in 
decibels is stored in YMDB. The three vectors YP, YM, and 
YMDB are all plotted with subroutine PLTTRF. The Versatec 
plot of the unit circle is generated with subroutine UNTCR. 
The 4912 plot of the unit circle and command array are 


generated with the subroutine UNTCIR. 


E. SUPERVISORY ROUTINES 

There are three supervisory routines used with the ASDF 
programs. These supervisory routines make the loading of 
ASDF transparent to the user. There is a substantial amount 
of file manipulation required to store and load _ the ASDF 
programs. A knowledge of the hardvare and software 
available for these purposes is not particularly useful to 
the user. For these reasons, the supervisory routines are 


discussed only briefly. Two of these routines are based on 
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CP/CMS and IBM EXEC environment; the third is written in the 
1BM=360 Job Control Language. 

The JCL program moves’ the LOUF = PEoechran from its 
permanent storage location on DISKO2 to the CP/CMS disks 
OSX001. The two EXEC routines generate an appropriate 
CP/CMS work space, and manipulate the ASDF text file into a 
usable form. 

The EXEC program ASDF191 resides on the user 191-P-disk. 
This routine automatically requests and formats ten 
cylinders of temporary disk space. The ASNF TEXT file is 
moved from the OSX901 disk to the newly created temporary 
disk space 192-T. Once in place, The ASDF TEXT file name is 
changed to ASDF EXEC to generate the next file with the 
correct file type. ASDF EXEC is then edited to remove the 
First set of data cards which hold the ASDF192 EXEC program. 
The 192-T-disk is released and the user is requested to type 
two commands to the system. The "login 192 P'" command makes 
the 192 disk the P-disk. The "EXEC ASDF192" command begins 
the execution of the second EXEC file. 

The general ASDF192 EXEC structure is depicted in figure 
3-6. ASDF192 EXEC is created by ASDF191 EXEC and resides on 
disk 192. In a manner’ similar to ASDF191, ASDF192 splits 
the ASDF TEXT file, which contains all of the component 
parts of the ASDF program, into a series of smaller files. 
Once all files have been separated, titles are altered 


appropriately, text files are loaded into core, and _ load 
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Figure 3-6 Structure Of ASDF192 EXEC 
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modules are generated. ASDF192 then reverts to an infinite 
loop, reading and implementing the desired user command. 
There are six user commands available to the user: POLZRO, 
RESPONSE, EXEC PUNCH, EXEC HRDSCPY, OPTIONS, and INTRO. The 
applicability of these commands is’ briefly discussed in 
chapter four. This brief description is also presented by 


executing the OPTIONS command once ASDF has been loaded. 
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IV. STRUCTURE OF ASDF COMMANDS 


The ASDF program generates the time and frequency 
responses of realizable digital filters of order ten or 
less. Chapter three describes the algorithms used for these 
computations. The procedure for executing the ASDF program 
generally consists of two distinct steps. The first step is 
to enter the poles and zeros of the filter to be analyzed. 
The second step is to perform the calculations required to 
generate and plot’ the desired filter responses. There are 
six commands’ which may be executed within the ASDF 
environment providing the user with the ability to control] 
the form of the output, to enter the poles and zeros of the 
filter, and to clarify most questions pertaining to’ the 
meaning and use of the zeneral commands. The six 
fundamental commands avallable to the user are: OPTIONS, 
meno, POLZRO, RESPONSE, EXEC PUNCH, and EXEC HRDSCPY. Only 
the POLZRO command has a_ subordinate command structure. All 
of these commands are issued by typing the command on the 
terminal keyboard in response to: *ASDF-RE. 

The OPTIONS command provides a listing of al] commands 
available to .the user. Each command is accompanied by a 
brief description of its intended use. 

The INTRO command writes on the screen six pages of 


information pertaining to the POLZRO command substructure. 
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INTRO includes a general description of the data input 
modes, and a description of several input situations with 
the appropriate executable commands. 

TRe POLZRO command allows the user to modify or create a 
filter in the z-plane. All adjustment of poles = and zeros 
must take place with this command. By issuing this command 
the user may add and delete poles and/or zeros as well as 
adjust the filter gain constant. The commands which are 
issued within the POLZRO environment are described in a 
later section of this chapter. 

The RESPONSE command displays the characteristics of the 
filter being created by the user. Five displays’ are 
presented: unit sample response, unit step response, phase 
of the transfer function, magnitude of the transfer 
function, and the magnitude of the transfer function’ in 
decibels. 

The EXEC PUNCH command allows the user to punch a card 
deck of the filter pole and zero locations on the offline 
punch. This punched output can be picked up in room 1[{-140 
of the computer center. The cards punched by this command 
may be subsequently reentered into the system so that this 
particular filter can be referenced or modified at a later 
date. 

When the user desires to enter a filter in this way, the 
punched cards are presented to one of the computer operators 


over the counter in room 1{-1490. The cards are punched in 
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the correct format for direct submission into the CP/CMS 
card reader. The operator subsequently reads this deck into 
the user's virtual card reader. As the user signs onto a 
terminal and gains access to CP/CMS, the computer includes 
an additional message in the standard login typeout showing 
that the user has a virtual reader file. This message is of 
the form: 
mee o esol RDR, NO Pai, NO PUN 

When the user wishes to access’ this old filter data he must 
execute the CP/CMS command OFFLINE REAN * *, The OFFLINE 
READ command will read the punched card filter data from the 
user's virtual card reader, and generate a file on the 
user's disk space entitled FILE FTOI1F901. Any files with 
this name present before the OFFLINE READ command is issued 
will be over-written by the new file. If the user does not 
wish to lose anolder FILE FTO1LFO001, he should use the 
CP/CMS ALTER command to change the name of the old file that 
he wishes to protect. Once the user has successfully loaded 
the ASDF program, the ASDF signature i.e., *ASDF-R8£&, is 
presented, The OFFLINE READ command may be issued any time 
that this signature appears. 

The EXEC HRDSCPY command is available only when the 
SUBMIT function of CP/CMS is functioning. When the SUBMIT 
command is not functional a message is sent to the’ user 
during the login procedure which states that SUBMIT is not 


available. the output from this command ts ultimately 
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available in the computer center, and consists of two parts. 
The first part of the output is the numerical listing of the 
values of the responses computed by the ASDF program. The 
EXEC HRDSCPY command is the only command which will allow 
the user access to these numerical values. The second part 
of the output is a series of Versatec plots, each similar to 
the screen presentations of the RESPONSE command, with one 
additional plot of the unit circle and pole/zero locations. 
When the EXEC HRDSCPY command is executed the ASDF program 
will ask several questions, using the responses to generate 
a JOB card. The response to the "ENTER DESIRED JOB NAME" is 
the title which appears on the final output. ASDF then 
submits the filter characteristics to the operating system 
for processing. The hard copy output will be found in the 


standard output boxes. 


A. SUMMARY OF DATA !NPUT MODES 

There are three modes through which the poles and zeros 
may be entered into the z-plane: 1) interactive graphics, 
2) rectangular, and 3) polar. Each of these modes’ has 
Particular attributes and limitations. 

1. Interactive Graphics (IAG) 

This mode was designed to enter poles and zeros 

rapidly, in cases where accuracy is not of prime importance. 
The interactive graphics mode is characterized by’ the 


appearance of the graphics cursor (the intersection of two 
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lines which trace across the terminal screen). The position 
of the cursor is controlled by the two thumb wheels located 
to the right of the interactive ASC!I! keyboard on the 
Tektronix 4012 terminal. The coverage of the interactive 
graphics cursor is limited to values of X and Y greater than 
-1.1 and less than +1.1. Zeros outside this area must be 
entered with either the polar or _ rectangular commands. 
Careful adjustment of the cursor can usually locate poles or 
zeros within .002 units of the desired location. Poles must 
be located within the unit circle. 
2. Rectangular 

This mode was designed to enter from. the terminal 
keyboard, poles and zeros which are known in” rectangular 
coordinates. It Is slower than the interactive graphics 
mode because there is ae significant amount of § screen 
rewriting required. The rectangular mode allows = much 
*sreater accuracy in the location of voles and zeros. Poles 
may not be entered on or outside the unit circle. Zeros must 
be located within ten units of the z-plane origin. 

3. Polar 

This mode was designed to enter from the terminal 
keyboard, poles”7 and zeros which are known in polar 
coordinates. This mode is also slower than the interactive 
graphics mode. t!{t allows poles and zeros to be entered to 
any permissible position in the z-plane. Poles must be 


entered within the unit circle. Zeros must be entered within 
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ten units of the z-plane origin. Each polar command has a 
description of the constraints on the data to he entered 
automatically written on the screen. This mode is 
particularily appropriate when constructing all-pass linear 


phase filters. 


B. THE POLZRO COMMAND SUBSTRUCTURE 

In an attempt to maximize execution sneed, a minimum of 
information is written on the screen, thus reducing the 
rewrite time after execution of each command. When program 
execution begins, the user is presented with the unit 
circle, and a listing of commands on the top left side of 
the screen. The commands are abbreviated in’ the following 
manner. The general command type or mode is-~ specified 
across the top of the array as INTERACTIVE, RECT, and POLAR 
(see figure 4-1). Each of these modes conforms to the above 
description. Down the left side of the array are_ the 
commands of the form CR8R, where C is the change to be made, 
A - for adda root, D- for delete a root. R&R is for the 
type of root to be acted upon by the command: RZ =o real 
zero, CZ - complex zero, RP - real pole, and CP - complex 
pole. As an example, the box in the center column alizned 
with the row titled DCZ stands for the command "Delete a 
Complex Zero - Rectangular Coordinates. There are two 
unique commands listed as NUC, and LST. The New Unit Circle 


(NUC) command simply redraws the unit circle and command 
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Figure 4-1 Basic Display Format 





array. This is appropriate when either too many” entries 
have been made to the terminal, and the command array is no 
longer clear, or when noise in the line from the computer 
causes a deformed array or unit circle to be drawn. The LST 
command provides the user with a listing of the current 
status of the filter being analyzed, in a self explanatory 
format. 

In all cases when the unit circle and command array just 
described appear on the screen, an audible tone or bell is 
heard. This bell indicates to the user that the cursor is 
to be positioned within one of the boxes in the command 
array. In order to execute a command, the cursor is 
positioned in the box associated with the desired command. 
A single keyboard character (any character but the "return" 
key is acceptable) is then typed. The cursor. then 
disappears. The system presents a carat ovdointing to. the 
command that the computer is about to execute; the flashing 
cursor blinks within the box associated with the command 
which is about to be executed. I{f the box so indicated is 
the command that was selected, the user types a "space", an 


matt 


(for execute) and then "carriage return". In cases 
where a rectangular or polar mode command is to be executed, 
the screen is then erased, and new information is printed. 
In cases where an interactive eraphics command is to. be 


executed, only the cursor reappears. This second appearance 


of the cursor is not accompanied by a bell or audible tone. 
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This indicates that the cursor is to be positioned within 
mer unit circle. Once so nositioned, a single keyboard 
character is typed (not "return'), and, as appropriate, the 
unit circle is redrawn to reflect the change in pole or zero 
locations. After this update of the z-plane, the circle and 
command array are redrawn with the cursor accompanied by the 
audible tone. 

There are only two exceptions to this scenario. First, 
if the flashing cursor does not flash in the correct box, 
the wrong command will be executed. If an incorrect box 
contains the flashing cursor, the user should not type the 
"space", "e'', sequence. Instead, he should type two other 
characters (e.g., "space", "space") and then "return." The 
cursor returns along with the audible tone, calling for a 
readjustment of the cursor, closer to the desired command 
box. The second exception occurs when the placement of al] 
of the poles and zeros is completed. In this case, after 
aligning the cursor with any box, the sequence "space", "x" 
(for exit), “return''’ should be typed. This will terminate 
the input, or modification, stage of the program and present 
the user with a series of options for the next’ stage of 


filter processing. At this point adisplay appears which 


explains the options available to the user. 
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C. SUMMARY OF POLZRO SUBCOMMANDS 
When the cursor appears along with an audible tone, the 


user is to position the cursor within the command array, and 


type a keyboard character (not "return''). When the flashing 
cursor appears within a box, he types ("'space", "tel! 
“return'™) to execute the specified command. "Snace'' and any 


other character but "x" will return the user to the command 


table for cursor readjustment. To terminate the input, or 
modification mode, the cursor is positioned over any box, 


and once the flashing cursor enters the box, the characters 


merece’, "x", “return” are typed. 
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V. LOADING THE ASDF PACKAGE 


The procedure for using the ASDF program consists’ of 
either six or seven steps depending on whether or not the 
ASDF package is resident on OSX901. Usually the program 
sections of ASDF are stored on the T2314 disk pack OSX901, 
Each Monday morning the OSX901 disks are scratched so that 
the ASDF program must be moved from the permanent storage 
location on DISK92, to the temporary storage location on 
OSX001. The user may simply sign onto a 4012 terminal, and 
execute three commands to load ASDF into his” virtual 
machine. Once this loading process has’ been successfully 
completed, the commands opreviously discussed are at’ the 


Mees disposal for filter analysis. 


A. LOGIN PROCEDURE 

With one minor exception, the standard losin procedures 
(specified in Section &.4.3 of the CP/CMS User's Manual) 
(Ref. 4] are applicable when entering the CP/CMS system to 
use ASDF. Due to the size of the ASDF programs” the user 
must request a 400k virtual machine when loging onto CP/CHS. 
The virtual machine is requested as follows: 

login XXXXTNN &O0k 

where XXXX is the user's computer center number, T is the 


type of user i.e. G for general user, and P_ for private 
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user, and NN is’ the two digit terminal number. The 4&00k 
request in the LOGIN is absolutely mandatory. Without this 
request the user will be able to load the ASDF program, 
however, when attempting to execute the RESPONSE command, an 
error: "INPUT COMMAND ERROR" will result. Should this error 
occur when the command is correctly spelled the user has no 
alternative but to logoff and log hack on to the’ system 
requesting the correct size virtual machine. The remainder 
of the login procedure is identical to the examples in the 


CP/CMS user's manual. 


pee tESt!NG OSX001 FOR ASDF PROGRAMS 
The second step in loading the ASDF programs is-7~ to 

determine whether or not a copy of the programs exists on 
OSX001. The ASDF programs are permanently stored on the 
operating system DISK02. The first time that the programs 
are used each week, the user must move them from DISK92 to 
the disk OSX901. To determine whether or not OSX901 has the 
ASDF package, issue the command: 

HOC each 1S) EXEC ert 
The system will respond: 

EXECUTION BEGINS... 

237 = OSX901 

ENTER DSNAME OR '/«! 
The correct user response is: 


FO718.asdf191 
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At this time one _ of two responses will appear. The 
particular response received is dependent on whether or not 
ASOF resides on OSX901. 

ee ASDF is on OSX001 

When the ASDF programs are located on OSX901, the 
system response is: 
POPENUoeeAnAUSTEN, EOF SIMULATED. 
4 BLOCKS 63 RECORDS READ 

This message indicates that File ASDE L911 has 
satisfactorially been transferred to the user's virtual 
machine and the user is ready to continue with’ the 
subsequent loading steps. 

2. ASDF is not on OSX901 

When the ASDF programs are not on OSX901, the system 
response will be: 
DATA SET NOT FOUND (0S/360 S213) 
ENTER DSNAME OR '/«! 

This response indicates that the user must move’ the ASDF 
programs to OSX001 before the loading procedure may be 
resumed, The simplest method for accomplishing this 
transfer is to generate a file on CP/CMS which exactly 
duplicates the Job Control Language program which follows. 

foe Zenerate the fite which transfers ASDF, logon to 
CP/CMS. Once the system's responses have been entered and 
the system is ready, issue the command: 


EDIT ASNF JCL 
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Then type in, or read ia the program which follows. 


Peeestandard JOB Card (User's Manual Sec. 3.3.2.1) 

// EXEC PGM=!IEBGENER 

J/SYSPRINT DD SYSOUT=A 

J/SYSIN DD DUMMY 

//SYSUT1 DD DSN=F0718.ASDF,VOL=SER=DISK02,UNIT=3330, 

// DCB=(RECFM=FB,LRECL=80, BLKSI ZE=3200) ,DISP=(OLD, KEEP) 
//SYSUT2 DD DSN=F0718.ASDF,VOL=SER=NSXON0L,UNIT=2314, 

if DCB=(RECFM=FB, LRECL=80,8LKSIZE=3200),D0ISP=(NEWN, KEEP), 
// SPACE=(TRK,(40,1),RLSE) 

// EXEC PGM=!EBGENER 

//SYSPRINT DD SYSOUT=A 

/J/SYSIN DD DUMMY 

//SYSUT1 DD DSN=F9718.ASDF191, VOL=SER=DISKN2,UNIT=3330, 
// DCB=(RECFM=FB, LRECL=80, BLKSI ZE=800) ,DISP=(OLD, KEEP) 
//SYSUT2 DD DSN=F9718.ASDF191, VOL=SER=O0SX001,UNIT=2314, 
// DCB=(RECFM=FB,LRECL=80,BLKSIZE=-800),0I1SP=NEW, KEEP), 
nf SPACE=(800,(3,1),RLSE) 

// 


Once this program has been correctly read in to’ the 
user's disk, it must be submitted to the operating system 
using the standard SUBMIT command: 

SUBMIT ASDF JCL Pl XXXX 
XXXX is the user's computer center number. If the SUBMIT 
Function is working the following messages will be returned: 
** CARDS XFERED TO 0098P a 
FNAME FTYPE Pl SENT FOR OS BATCH JOB 
After a short delay the system will send a second message: 
FROM 0998P33 : JOB SENT TO OS WITH 21 CARDS WRITTEN (Time) 


The user should then revert back to section B above which 


aA deck containing this program is available from 
Professor Kirk, Sp 302. 
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discusses the TOCP command. After giving the system a period 
of time to execute the JCL job which was just submitted, 
execution of the TOCP command will indicate that the ASDF 
files are present. If the operating system has not completed 
executing the transfer, the TOCP request will respond as in 


paragraph two above. 


B. LOADING ASDF 
Oncee ASDF resides on disk OSX901, the remaining required 
Steps are simple. The user issues the command: 
eee ao) Eola 
The system will type a series of lines onto the terminal. 
At the conclusion of these typewritten lines the user is 
prompted to issue the next two commands: 
LOGIN 192 P 
The system responds: 
Poe REPL AGEs |r 1491) 
The next required entry is: 
PREG wa see 192 
These two commands will cause a very long list of CP/CMS 
commands to be typed on the screen. Each of these CP/CMS 
commands is automatically executed by the operating system 
in the process of loading the ASDF programs. There are no 
furthrer commands required from the user to complete’ the 


loading process. 
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‘hen the entire ASDF program nackage ts loaded, the user 
will be presented with a title page and the tnstruction: 
NHIT SPACE AND RETURN TO CONTINUE." Once the user has typed 


Mepace' and “return” the display will write out a listing of 


the commands at his disposal. These commands have already 


been discussed in chapter four. 
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Vi. USING ASDF 


Each of the following sections constitutes an example 
for one of the three input modes. Each example uses only 
one input mode. In general, these modes may be combined 
during the generation of any particular filter. The mode 


selected should be the one most convenient to the user. 


PeeeeAAMPLE ONE - INTERACTIVE GRAPHICS 

This first example uses’ the ASDF program in_ the 
interactive graphics mode. The objective is to observe the 
responses of a digital filter which has two real zeros, one 
at the origin and another at +0.8, two complex poles located 
meme. cs t/~- j0.2), and the filter gain ts 1.000. The 
interactive graphics mode is used since the exact locations 


of the poles and zeros are not critical. For the purposes of 


this example, it is assumed that the user has’ already 
followed the login procedure described in the above 
Paragraphs, and has’ loaded the ASDF program. The user has 


already seen the title page, and the two pages of command 
descriptions. The screen is clear except for the top left 
corner which shows: 

AO Deion 


This group of letters is the signature for the ASDF program 
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and indicates that the program is ready to execute’ the 
user's next command. 

Since we must first enter the poles and zeros into the 
z-plane the appropriate command to type is: 

POLZRO 

The ASDF program will then draw the screen as shown in 
figure 6-1, and generate an audible tone. Initially the 
cursor will be located at some arbitrary position or perhaps 
may not appear on the screen at all. Movement of the thumb 
wheels moves the cursor position. Since we wish to execute 
in the interactive graphics mode, to add areal zero, we 
move the cursor (intersection of the two lines” tracing 
across the screen) into the box in the left column’ labeled 
ARZ, (Add a Real Zero) figure 6-2. We then type. any 
keyboard character except “return! > The trace cursor will 
disappear, and the flashing alphanumeric cursor will appear. 
The program will print a carat pointing to the box selected 
by the user. The flashing cursor will move into the box 
selected. 

The user then types the command "snace", "el", "return". 
The flashing cursor disappears, and the trace cursor 
reappears. The absence of an audible tone at this’ time 
indicates that the trace cursor is to be positioned within 


the unit circle. 


lL 7 : : 
The &012 terminal in the computer center requires a 
"control-s" after every response. 
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The trace cursor is adjusted to the nosition of the desired 
real zero, (the origin in this case), as shown in figure 
6-3. Any keyboard character except "return'' is then typed. 
The ASDF program onlots the zero at the origin, and returns 
the trace cursor with the audible tone. The audible tone 
reminds the user that the trace cursor is to be positioned 
in the command array. Since another real zero is to be 
entered, the cursor is positioned as in figure 6-1b again, 
and the same steps are followed. When the trace cursor 
reappears without the audible tone, position it within the 
unit circle as near as possible to the point (0.8 +/- j0.0) 
as shown in figure 6-4. (To enter the zero precisely at 0.8 
the rectangular or polar mode would be used.) The steps for 
entering the complex pole pair exactly parallel those for 
the previous entries, and are depicted in figures 6-5, and 
6-6. 

Once all the poles and zeros have been entered, the 
trace cursor may be aligned with any command box, and a 
character typed (not "return''). Once the alphanumeric 


oe 


cursor enters that box, the user types "space", 
"return''. The screen is cleared and the user is asked if he 


wishes to change the default filter gain. Since the default 


value is 1.0, the user responds "n", the ASDF program clears 
and responds with: *ASDF-RE 
The next step in the normal sequence is for the user to 


look at the filter characteristics. To display these 
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Trace Cursor Position Required To Place A Real Zero At The Origin 


Figure 6-3 





OL | oemeite ie SC MrINlDlS 


ocx 
ie 
rs 
= C 
mo 

_— 

Lib 
SO) 


SS eee ee 


aes I 


ae ee eae ee 


Pee Nae (OP a Cor 
ee Ge cee tec am) 
GE GCE GE SE (Ga) aa) aa, aaa | 


36 


Trace Cursor Position Required To Put A Real Zero At (+0.8) 


Figure 6-4 
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characteristics the command typed is: 
RES POHSE 
To set a hard copy of the filter responses the command: 
EX EGmEROSCPY 
is typed. Figures 6-7a, b, c,d, e, f, and 2g show. the 


pertinent portions of the output. 


Bec X AMPLE TWO = RECTANGULAR MODE 

This second filter example is taken from (Ref. 5], page 
222. The filter is a fourth-order low vass Chebyshev. The 
derivation of the transfer function is described in’ the 
text. The given transfer function is factored to give the 
desired pole/zero locations. The filter has two complex 
zeros at (-1.0 +/- j59.9) and two complex pole pairs, one 
pair at (.7498 +/- 5.5348), and one pair at (./774 +/- 
Geel20). The filter gain is given as 0.001836. 

In @ manner exactly paralleling the previous example, 
the user selects the command ACZRCT (which stands for Add a 
Complex Zero - Rectangular Coordinates), by positioning the 
trace cursor in the appropriate box, as in figure 6-8a. 
Once a keyboard character is typed, and the alphanumeric 


Vee and “return 


cursor moves into the box, the "space", 
keys are typed. The ASDF program will erase the screen and 
draw the display shown in fisure §-8b. The user then tyoes 


in the rectangular coordinates of either complex zero in 


the first complex zero pair, figure 6-8c. The system wild 
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ask whether or not the root location has been correctly 


entered. The user responds 'y' for yes, 'n' for no. If the 
root has been incorrectly entered, an 'n' response will 
rewrite the screen for another attempt. Upon receiving a 
'y' response ASDF rewrites the unit circle and command 
array, and plots the new root. The second complex zero is 
entered exactly as the first. 

The complex poles are entered as were the zeros. The 
trace cursor is first placed in the box ACPRCT. (Add a 
Complex Pole - Rectangular Coordinates), as in figure 6-9. 
Figure 6-10 shows the unit circle after all poles and zeros 
have been entered. The POLZRO command is terminated by 
aligning the trace cursor with any box, typing a character, 
waiting for the flashing alphanumeric cursor to enter the 
box, and typing "space", "x". The system responds by asking 
about the filter sain. In this example the filter gain is to 
be changed to .001836, so the correct response is 'y'. The 
new filter sain, 9.001836, is typed to terminate the POLZRO 
command. 

To view the filter responses, the user enters’ the 
command: 

RESPONSE 
when the *ASDF-RE signature appears. The hard copy 


responses of the filter are shown in figures 6-lla through 


6-lle. 
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Unit Circle After All Poles and Zeros Have Been Returned-Example Two 


Figure 6-10 
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Bee GXAMPLE THREE - POLAR MODE 

This example uses the polar mode to Zenerate an all pass 
Filter. The filter complex pole pairs are to be located at 
(09.2 at +/- 4&5 degrees), and (9.5 at +/- 30 degrees). The 
filter zeros are to be at (5.0 at +/- 4&5 degrees), and (2.0 
at +/- 30 degrees). The filter zain is to be 90.01. The 
first step is to change the specifications in degrees to 
radians. Thirty degrees equals .5235 radians. Forty five 
degrees is .7854 radians. The entry of ooles and zeros 
exactly parallels the above examples, however, the desired 
comands are ACPPLR&, (Add a Complex Pole - Polar 
Coordinates), and ACZPLR, (Add a Complex Zero - Polar 
Coordinates). Figure 6-12a shows the trace cursor position 
mom ACPPLR. Figure 6-12b shows the trace cursor position 
for ACZPLR. The angular part of any pole or Zero entered in 
the polar mode must be in radians. Figure 6-12c shows the 
entry of the first complex zero pair. Figure 6-12d shows 
the entry of the first complex pole pair. Figure 6-12e 
shows the change in filter zatin from 1.00 to 0.01. Figures 
6-13a through 6-13g show the hard copy responses’ for this 
Filter. !t should be noted that zeros which plot outside of 
the unit circle are indicated by stars. All such zeros are 
Plotted at a radius of 1.1 along the same angular component 
as the actual zero location. Zeros are plotted this way for 
two reasons: first to make the output display consistent for 


all relevant filters, and second to emphasize zeros which 
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Figure 6-12b Trace Cursor Position For The ACZPLR Command 
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are close to, but outside, the unit circle. The importance 
of distinguishing between zeros very near the unit circle, 
either interior or exterior, rests on the extreme 
sensitivity of the phase of the frequency response as a zero 
moves from the interior, onto, and to the exterior of the 


mae circle. 


a3 





YVIT. “ADDITIOHAL EXAMPLES 


The purpose of this chapter is threefold. TMS acieY 
objective is to provide the user with a series of plots 
demonstrating the general form of those types of filters not 
considered in chapter. six. Additionally, this chapter 
denonstrates the responses generated by poles located at a 
series of positions within the unit circle. It further 
provides examples for user comparison when learning to 
interact with the z-plane. 

The filters considered in the examples of chapter six 
include a fourth-order low-pass filter and a second-order 
all-pass filter. The sequence of figures 7-la through 7-1¢ 
shows the general pole zero placements for ae high-pass 
Filter, and associated responses. Similarly, fiture 7-2a 
throush 7-22 show the responses for a band-pass filter. 
Figures 7-3a through 7-32 show a band-reject or notch 
filter. 

The secondary purpose of the chapter is to show the user 
the general forms of time and frequency responses as the 
filter poles move’ in the z-plane. The sequence of figures 
7-4 through 7-11 shows the responses for a_ real pole at 
various locations on the real axis. Figures 7-12 through 
7-19 similarly show the responses of a complex pole pair 


movine across the 2-plane. Frequency responses are not 
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shown for poles on, or outside of, the unit circle since the 
frequency response is not defined for these cases. Figures 
7-17, 7-29, and 7-21 show the change in responses for a Pair 
of complex poles which move away from the real axis. 

The final Intent of the chapter is to provide the user 
with additional figures for comnarison should he _ need 
further practive in manipulating poles and zeros in_ the 
z-plane. It is itmnortant to note that the ASNDF nrosram wil) 
not allow the user to position poles on or outside the unit 
circle. This limitation results in several plots which can 


not be reproduced by the user. 
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Vt weenie toner —ASDE OUTPUT 


When the user initially confronts the ASDF program it 
is expected that he will attempt to duplicate some of the 
examples presented earlier’ in this thesis or similar 
examples found fn the literature. These tyoes of 
applications of ASDF are straightforward and the’ results 
zenerated can be compared with published results. 
Interpretation of the output plots is seldom confusing. 

As the user becomes more familiar with the use of the 
ASDF program and the positioning of the poles and zeros in 
the z-plane he will undoubtedly attempt to generate filters 
with sharper transition regions, and greater attenuation in 
the stop band(s). In general these filter improvements are 
accomplished by increasing the order of the system = and 
positioning the poles closer to the unit circle. As the 
order of the system increases, and as the poles of the 
filter become closer to the unit circle, the filter becomes 
increasingly subject to the difficulties associated with 
the finite precision arithmetic used in the filter 
implementation. 

A tacit assumption made very early inmost digital 
signal processing courses (and often quickly forgotten) is 
that infinite precision arithmetic is available in 


performing the calculations. This assumption breaks down 
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when digital filters are realized in hardware or simulated! 
ln general there are three sources of error which may affect 
the performance of digital filters. Due to the finite 
precision representation of the coefficients, the values of 
the implemented filter coefficients are only approximately 
equal to the desired values. Second, errors may be generated 
in the arithmetic operations. In general, the product of 
two n-bit numbers results ina 2n-bit result. An error 
occurs when this 2n-bit number is converted to n-bits for 
succeeding computations. The third type of error 
encountered is the analogf to digital conversion required for 
processing analog signals with dirital filters. Only the 
first two of these errors is of concern in the ASDF program. 

While an in depth analysis of finite precision effects 
is beyond the scone of this work, let us consider the filter 
which has eight zeros located at the origin, and eight poles 
located at -0.9. Entering these locations into the z-plane 
using the ASDF program generates the filter responses’ in 
fizures 6-la, b, c. A superficial glance at these responses 
indicates that this filter is unstable. The filter must be 
stable, however, since all of the poles lie within the unit 
circle. 

There are several phenomena involved in this example. 
The first is the dynamic range limitation of the system on 
which the filter is implemented. The I[BM-360 is limited to 


computing with numbers between Ga and ice ~. AYnce the 
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larger of these limits is exceeded, the computations are 
terminated and it iis not clear whether’ the designed filter 
is actually unstable or has merely exceeded the dynamic 
range limitations of the machine. The dedicated user may 
decide to modify the programs presented in the appendicies 
to investigate the time responses of such systems. Clearly 
the form of the time responses will not’ be significantly 
altered if a modified unit sample response 

Mr ed, (0, 05 « -« 
and a modified step response 


iis: OM en Ole) x 10) 


Ie 
, 
is used. 
if these modifications are nerformed, two other 
phenomena become apparent, each relating to the finite 
precision implementation of the digital filter. The first is 
that, in general, digital filters have deadbands, the second 
that the filter output may degenerate to a limit cycle. The 
original works by Blackman [Ref. 1] . demonstrate that 
digital filters exhibit both of these effects. Extensive 
covgerage of these phenomena is not appropriate here, 
however, a dead-band in a digital filter can be 
simplistically described as a span of filter input values 
for which the filter output does not change. The limit 
cycle effect is evident when a digital filter responds to a 


steady state input. For constant inputs, the filter output 


May run through a cycle of values. The pronensity for 
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experiencing limit cycles is obvious in the calculation of 
the step response. A close look at the unit sample response 
shows that after a sufficiently larze number of input 
values, this unit sample sequence becomes in reaching 
steady state, a constant. 

When observing the output response of the examnie 
filter, there is still another possibility for explaining 
this output. This explanation revolves around the fact that 
the poles in the implemented filter are not at the locations 
expected. In fact, the poles may actually move outside the 
unit circle even though they were entered within. In 
general, there are only a finite number of pole locations 
possible with finite word lensth coefficient 
representations. The effect of this limitation is best 
explained using another example. Consider a filter with 
three zeros located at 9.5, and 0.2 +/- J0.5, and four poles 
located at -0.999, -0.6, and -0.69 +/- J0.69. The infinite 
precision transfer function for this filter is expressed in 


equation 8.1. 


Peas 3 Fe eee. TAS C2) 


3 + 3.75322 2° + 2.3497398 Z + .57074868 





me 2,979 2 


Now consider the finite precision case where each of the 


coefficients must be represented by a six bit number. The 


pole at -9.6 is expressed in binary (exclusive of sign) as 
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meoeloOTTO9. . . and is truncated to six bits as .100110. 
When all of the coefficients are so represented, the 
transfer function is expressed by equation 8.2. Clearly the 
poles and zeros of this second transfer ftunction are not the 


7° = .890625 2° +.474275 Z - .140625 (8.2) 


H(z) = ES = 7. 
Peerage t+ 5.7502 ° + 2.5125 2 + 1.5625 


same in equation 8.1. In general the differences will change 
all of the filter responses. In fact for this particular 
case the new poles of the system are at -9.79522 +t/- 
J0.16677, and -9.76353 +/- JO ,o9cotoyweenis ~§ filter --ts 
unstable! 

The reason for the shift in pole’ locations Is clearly 
portrayed by the filter whose characteristic equation Is 
equation 8.3. 

(Z + .995)% = 9 (8.3) 
This filter has eight poles at -0.995. Note what happens to 
the filter when an error is introduced by truncating the 
infinite precision coefficient values to thirty-two bits. 
For simplicity restrict the discussion to the ramifications 
of an error generated in the value of the constant (.995)° 
This error can be as large as 2 xX Wome Wrehe ents Error 


explicitly [ncluded, the new characteristic equation will 


he: 
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The roots of equation (38.4) are complex, and lle ona circle 
of radius 0.06 around the point -01.995 #/- J 0.0. Pier 
8-2 shows that three of these roots lie outside the unit 
circle. Thus with a finite precision implementation of the 
filter (even with 32 bits of precision), it is possible for 
the pole locations to move so far that the filter actually 
becomes unstable. 

The results of this section as well as four’ types of 
ASDF outputs which have been observed, are summarized in 
figure 8-3. Trace A shows the unit sample response of a 
stable filter. Trace B shows the response of a filter which 
exceeds the dynamic range limitations of the implementation 
system. Trace 81 is for an unstable filter and B2 is for a 
stable filter caught in a limit cycle. These two traces can 
be observed or explored only by modifying the unit sample 
Sequence. Trace C shows a stable filter which falls into a 
limit cycle below the dynamic range limit of the machine 
upon which the implementation has been performed. Note that 
the performance of a filter above the dynamic range 
limitation of the 18M 369 (Cor any other machine) cannot be 


observed, unless the ASDF program is modified. 


Z09 





eToOuT9 1Tun eu, eptsizno seuy SeTog xeTdwog esau, Jo eeuuy 





eToutg 4tTuy 


2-8 eUNnsTy 


Za 





Sosuodssy oetTdwes tug peAdesqg eTqtssog jo Aaewung E€-g sUunstTy 


od 


td 


LIN 
AIONVa 
JIWVNAG 


AISNOdSda 
ATdWVS LINN 
SO elileliy 2 


Za 





IX. SUMMARY AND CONCLUSIONS 


The Advanced Simulation of Digital Filters has been 
implemented on the 18M 360, and is intended for the use by 
persons beginning their study of digital signal processing. 
The programs included within the package have been designed 
to make the user's interaction with the z-plane as simple as 
possible. While the programs themselves are available for 
study, all of the computations are invisible to the user, in 
an attempt to develop the user's intuition of the discrete 
domain without overwhelming him with arithmetic 
computations. 

The ASDF package provides the user with an interactive 
method by which the filter nole and zero locations’ can be 
manipulated and the effects observed of such movement. ASDF 
can simulate filters of up to tenth order, and is punctuated 
with numerous comments and instructions which are to assist 
the user. The interface with the z-plane is available in a 
rapid interactive e¢raphics mode, as well as in the slower 
polar and rectangular modes, each providing gfreater position 
location accuracy. The filter responses generated by ASDF 
are presented asa series of displays directly on the 
Tektronix 4012 terminal. Versatec plot capability has been 
integrated into the system to provide high quality copies of 


the system responses. 
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mm CONCLUSIONS 

There are several points which ere worth considering for 
persons intending to nursue this or other areas of 
interactive systems development with the graphics terminals. 
While comments may apnear obvious they are stated to 
minimize the difficulties which can be encountered. Sip Sten we 
serfous attempt should be made to define the entire data 
structure before beginning the programming stage of the 
work. All phases of the effort must be planned before 
commencing with the actual programming to insure’ that the 
data structure selected is sufficiently teneral, and 
sufficiently flexible to endure the additions and changes 
which inevitably occur when initial objectives are 
overlooked or changed. Second, it is imperative to master 
all cf the software tools available and anplicable to the 
project. There are features of the PLOT 10 software which 
exist and were not used in the ASDF program. This 
shortcoming is due to insufficient a priori knowledge. For 
example, the writing of text oortions of the screen 
presentations should have been much’ simpler. Finally, it is 
important to research all of the attributes of the devices 
used for the final problem solution, and direct the solution 


to take advantage of all the system capabilities. 
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B. AREAS FOR FUTURE WORK 

There are several areas in which this program could be 
improved or extended. The possibility exists ior 
transferring the ASDF programs to the Tektronix 4082 
computing system. This move would avoid many of the 
shortcommings associated with the current IBM 3650 system. 
The user would be relatively immune to system crashes, and 
the constraints placed tupon him by concerning the 18M 
system's being operational and/or available. There are 
several disadvantages associated with this move. With the 
advent of the new computer system the majority of the 
current problems concerning system reliability, and 
constraints on available access hours should be eliminated. 
Furthermore, transferring these programs to the 4082 system 
will require significant rewriting to make execution 
possible on asystem which has amuch smaller virtual 
machine. Furthermore to fully benefit from the 4082 system, 
the interfacing routines should be rewritten to. take 
advantage of the 4082 refresh capability. It seems that all 
of the effort involved for this change is not worth ene 
expected improvements. 

The ASDF program is now useful for analyzing digital 
filters. lt seems likely that this capability could be 
expanded into a method for computer-aided design of digital 
filters. There are standard transformations which allow the 


filter designer to map analog filter requirements into the 
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appropriate digital requirements, and as well there = are 
transformations which make possible the computation of 
digital filter transfer functions from known analog transfer 
functions. There are also aevariety of optimization 
techniques referenced in the literature which could be 
implemented. 

The other area of obvious expansion for the ASDF program 
is to include’ the capability for oerforming spectral 
analysis. lt would be worthwhile (though perhaps not 
particularly involved) to include in the system the ability 
to calculate Discrete Fourier Transforms, and = graphically 
display the results. Additionally, it would be possible to 
demonstrate the correspondence between performing 
convolution in the time domain directly and overforming the 
convolution fn the frequency domain by multiplying the 
DFT's. The relation between the circular and linear 
convolution could also be graphically portrayed. 

As a final precautionary note, it is not pnossible to 
overestimate the length of time required to scenerate well 
documented, correct, useful orograms on an interactive 
basis. Changes in system parameters, and the limitations in 
System software can make the process of simplifying the 
user's interactive responsibilities a very difficult process 
for anyone not thoroughly acquainted with all available 
resources. The myriad of possible incorrect user responses 


which must be anticipated make the problem extensive. 
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merenoixX Az: PROBABLE SOLUTIONS FOR DIFFICULTIES 


SYMPTOM: RESPONSE command returns an error message upon 
execution. 

SOLUTION: Check the virtual machine size with the command 
"CP QF" to insure a virtual machine stze of 4N0K. If the 
user has only been allocated 256, logoff and log back on as 
directed in chapter five. Load ASDF as described in chapter 


six. 


SYMPTOM: User has itnadvertently entered the contro] 
program, ie., CP. 
SOLUTION: moe ine) and "return," Normal orogram 


execution should resume, 


SYMPTOM: The carat appears with the blinking cursor, but it 
is not located in the command array. 
SOLUTION: Move the trace cursor into the command array, and 


type a "space." Normal operation should resume. 


SYMPTOM: No terminal activity. 

SOLUTION: First, check to see whether or not the system has 
crashed. If the system has crashed, logon as in chapter 
five, and load ASDF as in chapter six. If the system has not 


crashed, move the cursor to the center of the screen and 
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Eyope "space" and "return.' tf normal operation has not 
resumed hit the "break" key. When the CP prompt appears, 


type "b."' Normal operation should resume. 


SYMPTOM: An error is detected in the execution of ASDF191 
EXEC indicating that no temporary disk space is available. 
SOLUTION: Check with the consultant or a systems 


programmer. 


SYMPTOM: The user can not escape from a delete pole or 

delete zero command in the interactive graphics mode (IAG). 

SOLUTION: Implicit in each request to delete from. the 

z-plane is a request, validating the change to be made. In 

the delete IAG mode the trace cursor returns after 

illustrating the root to be deleted. Type "y" if the correct 
1b at 


root was identified. Type n” if the wrong root was 


identified, and attempt the procedure again. 


SYMPTOM: The user did not enter the correct amount of data 
and the program terminated on a 217 error. 

SOLUTION: Type "b" and "return.' The program = should 
restart with the ASDF signature. At this point reenter the 
command desired. If executing the POLZ20 command it is 
advisable to execute the list (LST) command to check on the 


eurrent filter status. 
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APPENDIX B: IMPORTANT SUBROUTINES AND VARIABLES 
The subroutines and variables described in this appendix 


are applicable to all programs within the ASDF package. 


A. SUBROUTINE NAMES 
The subroutines within the ASDF package are divided into 
three major catagories: 1) Plot-10 software, 2) Versatec 
software, and 3) routines generated explicitely for the ASDF 
programs. 
1. Plot-10 Software Subroutines 
The Plot-10 subroutines are all adaquately described 


in reference 10, and are listed below. 


ERASE MOVABS ANMODE 
RECOVR HOME LINEF 
NEWLIN SHE NDO VIWENDO 
POLTRN MOVEA DRAWSA 
LINTRN DRAWR DRVWREL 
MOVREL VCURSR Beeb 
POINTA SGU SS INIT 
TA8HOR AB we SEERIC 
FIN PARLE. DRAWA 
DRYWABS 
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2. Versatec Subroutines 
The Versatec subroutines are listed below, and are 


explained in detail in reference 7. 


EO PLOT WINDOW 
SCALE AXIS SYMBOL 
LINE 


3. ASDF Subroutines 

Each of the following subroutines have been 
generated for the ASOF programs. 
CRRMMM - Add and delete poles and zeros from the z-plane. 
The change "C"' can be either "A" for add or "D"' for delete. 
maiemroots "RR" can be "RZ" for real zero, "CZ" for complex 
zero, "RP'' for real pole, or "CP" for complex pole. Mode 
"MMM" can be "IAG" for interactive graphics, "POL" for 
memes, or "RCT" for rectangular. 
FLTR2 - Prints a listing of the pole and zero locations of 
the filter being considered. 
ASUBK ~ Generates the ay coefficients required to evaluate 
the transfer function. 
BSUBR - Generates the b, coefficients required to evaluate 
the transfer function. 
ADJUST - Insures that the correct number of zeros located at 
the origin are included in the coefficients ay and b,.. 
RSPNSE - Computes the time responses of the filter. 
FRQNCY - Computes the phase and magnitude of the BROAe Fair 


function for the prescribed filter. 
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PRPPLT - Scales the ordinate data prior to plotting the 
responses. | 

PLTIMP - Plots the unit sample response. 

PLTSTP - Plots the unit step response. 

PLTTRF - Plots the phase of the transfer function as well as 
the magnitude in linear and decibel format. 

UNTCR - Plots the unit circle with the appropriate poles and 
zeros. 

START - Initializes all arrays and required variables. 
FINISH - Checks to insure that the constructed filter is 
causal, and generates a file to save the pole and zero 
locations. 

ANGLE - Draws the angle sign for polar commands. 

ERROR = Controls orofgram flow and error messazes for user 
generated errors. 

LST - Prints a listing of the current pole zero location. 
PLTTBL - Plots the poles and zeros on the unit circle. 
CHCKSZ - Insures that roots are available to be deleted and 
that the maximum system order will not be exceeded. 

UPDATE _ Makes changes in the POLZRO array to corresoond to 
user changes in the z-plane. 

LOCATE - Searches the POLZRO table for the root closest to 
that which the user wishes to delete. 

STRRTS - Puts poles and zeros into the POLZRO array. 


PLOTRT =- Plots a single pole or zero on the unit circle. 
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B. IMPORTANT VARIABLES 

The variables described below pertain to all of the 
subroutines and programs within the ASDF voackage. 
IROOTS(!) - Keeps track of the number of poles and zeros in 
the system under consideration. !ROOTS(1) is the number of 
real zeros. IROOTS(2) is the number of complex zero pairs. 
IROOTS(3) is the number of real poles. {R9O0TS( 4) Is the 
number of complex pole pairs. 
POLZRO(I,J) - The locations of the zeros are stored for 
values of !| from one to ten. The pole locations are stored 
for values of | from 11 through 20. The real part of the 
rectangular representation of a root location is stored with 
J equal to one. The imaginary part of the rectangular 
represention for a root location is stored with J equal to 
two. The magnitude and angle of the root's polar 
representation are stored for J values of three = and four, 
respectively. POLZRO(!,5) holds the root type. 
IERROR - Keeps track of the generation of an error. Once 
the error has been identified the value is set to one. When 
the error has been expnplained to the user the value is reset 
to zero. 
AC!) and BCI) - Store the values of the ay . 
coefficients. 
AAC 1) and BBCI) =- Store the double precision values of the 


ay and oe coefficients. 
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YICl) - Stores the values of the unit sample response to be 
plotted. 

fee, «stores the values of the unit step response to be 
plotted. 

YP(!) - Stores the values of the phase of the transfer 
function to be plotted. 

YM( I) - Stores the magnitudes of the transfer function to be 
plotted. 

YMDB(I) - Stores the transfer function magnitudes’ in 
decibels to be plotted. 

EN(I) - Stores the time values against which the dependent 
variables are plotted. 

IYIPTS - Holds the number of unit sample points to be 
plotted. 

IYSPTS - Holds the number of unit step response points to be 
plotted. 

NDEN - Is the order of the denominator plus one. 

NNUM - Is the order of the numerator plus one. 

NORDER - Its the order of the system under consideration. 
This can be either of the numberator or denominator. 

LAST - Holds the number of points to be plotted. 

FLTRGN - Is the value that the user assigns to the filter 
Zain. 

ZEROS(1) - Holds the complex coefficients which will be made 


real to generate the b,. coefficients. 
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POLES( I!) - Holds the complex coefficients which will be made 
real to generate the a. coefficients. 

The remainder of the variables are either obvious, or 
used merely for control of the program flow, ie. do loon 


indicies, etc. 
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Meet nk C. STRUCTURE FOR ASDF TEXT FILE 


The file ASDF TEXT is a composite file made up of all of 
the required programs and EXEC files for running the ASDF 
program package. These files have heen consolidated to 
minimize the time required for transfer from DISKN2_ to 


0SX001. The structure of the comnosite file is listed below: 


BREAKOO1 

EXEC file ASDF192 

BREAKNO2 

TEXT file for POLZRO command 

BREAK903 

TEXT file for RESPONSE command 

BREAKIO4 

Information printed out by the INTRO command 
BREAKI05 

meeecl JCL - used for EXEC HRNSCPY 

BREAKI06 

VTEC2 JCL - used for EXEC HRDSCPY 

BREAKI07 

Lables for RESPONSE plots 

BREAKI08 

EXEC file to execute the EXEC PUNCH command 
BREAK909 


224 





Initial file FTOLFO01 with no poles or zeros 
BREAKON10 

TEXT file for generating JO8B cards 
BREAK)11 

TEXT file for the INTRO command 
BREAKI12 

Information printed by OPTIONS command 
BREAK913 

Information printed on the thtle page 
BREAKON14 

EXEC file to execute HRDSCPY command 
BREAKI15 

TEXT file for the OPTIONS command 
BREAKO16 

TEXT file for the software page clear 


BREAK)17 





PPPENDIX D. JCL FOR EXEC HRDSCPY COMMAND 


The Job Control Language required for the execution of 
the EXEC HRDSCPY command is provided to sive insight into 
the file manipulation utilized by the ASDF package. Should 
the permanent locations of any of the portions of the ASDF 


package change this section of JCL must be changed also. 


GE for VTEC! 


//GO EXEC PGM=THEPLT,REGION=180K 

//STEPLIB DD UNI T=3330, VOL=SER=NISKN3,DISP=SHR,DSN=S1502.ASDF 
//FTO5FO001 DD DDNAME=SYSIN 

//FTO6FN01 DD SYSOUT=A 

//FT15F001 DD DDNAME=PLOTPARM 

//PLOTLOG DD SYSOUT=A 

Me ECIR2 DD DSN=%&VECTR2,N1SP=(, PASS) ,SPACE=(CYL,(5,5)),UNITSSYSDA 
//VECTR1 DD DSN=4&VECTR1,DISP=(, PASS), SPACE=(TRK,(1,1)),UNIT=SYSDA 
//GO.SYSIN DD * 


JCL for VTEC2 


//PLOT EXEC PGM=IEVMAPP, COND=(4,LT) 
//FT15F001 DD DDNAME=PLOTPARM 

MeerePLiB DDO DSN=SYS1.VTECPLOT,DISP=SHR 
//PLOTLOG DD SYSOUT=A 

//NECTR1 ODD ODISP=(OLD,DELETE),DSN=%&VECTRI 
//VECTR2 ODD ODISP=(OLD,DELETE),DSN=&8&VECTR2 
MmeorovECTR DD SYSOUT#=(A,,5555) 

//VECTTAPE DD DUMMY 

u/ 
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oo ey eae Once, DECK FOR ASDPF1L91L EXEC 
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